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Abstract 



£Sj . In this work we discuss the properties of the solutions of the Dirac equation 

in presence of an uniform background magnetic field. In particular we focus on 
the nature of the solutions, their ortho- normality properties and how these solu- 
fN) ■ tions depend on the choice of the vector potential giving rise to the magnetic field. 

We explicitly calculate the spin-sum of the solutions and using it we calculate the 
propagator of the electron in presence of an uniform background magnetic field. 

O ' 

1 Introduction 

Calculations of elementary particle decays and scattering cross-sections in presence of 
a background magnetic field are commonly found in literature [H El EH El IS]- These 
calculations became more important after it was understood that the neutron star cores 
can sustain magnetic fields of the order of 10 13 G or more. These realistic fields may be 
very complicated in their structure but for simplicity many of the times we assume these 
fields to be uniform. The advantage of an uniform magnetic field is that in presence of 
this field the Dirac equation can be exactly solved. Once the the Dirac equation is exactly 
solved then we can proceed to quantize those solutions and calculate elementary particle 
decays and scattering cross-sections in presence of the background magnetic field. In this 
article we will solve the Dirac equation in a background magnetic field and discuss about 
the nature of the solutions. We will quantize the fermionic theory in presence of a magnetic 
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field and will end with the calculation of the electron propagator. As the calculations rely 
heavily on the choice of the background gauge field giving rise to the magnetic field we 
will discuss about the gauge dependence of the various quantities calculated in this article 
and infer about the gauge invariance of physical quantities as scattering cross-sections and 
decay rates calculated in presence of a magnetic field. 



2 Charged fermion in presence of a magnetic field 

In this article we will assume that the uniform classical background magnetic field is 
along the z-direction of the co-ordinate axis. The background gauge fields giving rise to 
a magnetic field along the ^-direction, of magnitude B, can be fixed in many ways: 

A B = AI = A B = J Al = -yB + b. (1) 



or 



or 



A B = A x B = A z B = 0, A B = xB + c. (2) 



A° B = A B = 0, Al= l -xB + d 1 A B = - l -yB + gi (3) 

where b, c, d and g are constants. Here A B designates that the gauge field is a classical 
background field and not a quantized dynamical field. In the above equations x, y are just 
coordinates and not 4-vectors. In this article we will assume that the gauge configuration 
as given in Eq. (pQ) with 6 = 0. More over in this article we will be employing the 
Dirac-Pauli representation of the Dirac matrices. 

2.1 The solution of the Dirac equation 

The Dirac equation for a particle of mass m and charge eQ, in presence of a magnetic 
field is given by: 

<^ = W B ^, (4) 

where the Dirac Hamiltonian in presence of a magnetic field is given by: 

n B = d-U + Pm. (5) 

Here n M is the kinematic momentum of the charged fermion. In our convention, e is the 
positive unit of charge, taken as usual to be equal to the proton charge. From Eq. (jJJ) we 
can infer that for the stationary states, we can write: 

iP = e-^(^y (6) 
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where 4> and x are 2-component objects. With this notation, we can write Eq. (j3J) as: 



(E — m)4> 
(E + m)x 



Eliminating x-> we obtain 



(E 2 - m 2 



a- • (-iV - eQA) X , 
a ■ (-iV - eQA)<f>. 



-zV - eQA) 



With our choice of the vector potential, Eq. (J9j) reduces to the form 



2 m 2 ) 



V 2 + (eQSfy 2 - eQB(2iy^- + 03 



(7) 
(8) 



(9) 



(10) 



Here 03 is the diagonal Pauli matrix. Noticing that the co-ordinates x and z do not 
appear in the equation except through the derivatives, we can write the solutions as 



e^f(y) 



:n) 



where f(y) is a 2-component matrix which depends only on the ^/-coordinate, and possibly 
some momentum components, as we will see shortly. We have also introduced the notation 
X for the spatial co-ordinates (in order to distinguish it from x, which is one of the 
components of X), and X^ for the vector X with its y-component set equal to zero. In 
other words, p ■ X^ = p^x +p z z, where p x and p z denote the eigenvalues of momentum in 
the x and z directions!!! 

There will be two independent solutions for f(y), which can be taken, without any 
loss of generality, to be the eigenstates of 03 with eigenvalues s = ±1. This means that 
we choose the two independent solutions in the form 



U(y) 



F+(y) 




f-(y) 





F-(y) 



Since a 3 f s = sf s , the differential equations satisfied by F s is 

(/ ~ ''" - (eQBy + p x ) 2 F s + (E 2 - m 2 - p\ + eQEs)F s = , 



dy 2 



(12) 



(13) 



which is obtained from Eq. (flOl) . The solution is obtained by using the dimensionless 
variable 



Z = Je\Q\B y 



Pa 



which transforms Eq. fTTSl) to the form 



eQB 



F„ = Q 



(14) 



(15) 



1 lt is to be understood that whenever we write the spatial component of any vector with a lettered 
subscript, it would imply the corresponding contravariant component of the relevant 4-vector. 
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where 

E 2 -m 2 -p 2 z + eQBs 
tts ~ e\Q\B 



(16) 



This is a special form of Hermite's equation, and the solutions exist provided a s = 2v + 1 
for v — 0, 1, 2, • • -. This provides the energy eigenvalues 

E 2 = m 2 + p\ + (2i/ + l)e\Q\B - eQBs , (17) 

and the solutions for F s are 

N u e-?/ 2 H u (0 = I»(0, (18) 

where H v are Hermite polynomials of order is, and N u are normalizations which we take 
to be 

1/2 



'e\Q\B 
!2V^ 

With our choice, the functions l v satisfy the completeness relation 



"•=lfe=l • < 19 > 



Y. ««)/,«„) = sje\Q\B 5{i - {,) = % - y.) , (20) 

where is obtained by replacing |/ by in Eq. f|T4l) . 

So far, Q was arbitrary. We now specialize to the case of electrons, for which Q — — 1. 
The solutions are then conveniently classified by the energy eigenvalues 

K = m 2 + p 2 z + 2neB , (21) 

which is the relativistic form of Landau energy levels. The solutions are two fold degener- 
ate in general: for s — 1, v — n — 1 and for s = —1, v = n. In the case of n = 0, from Eq. 
f|T7|) we see that for Q = —1, v = — 1(1 + s), and as v cannot be negative s = —1. Thus 
the n = state is not degenerate. The solutions can have positive or negative energies. 
We will denote the positive square root of the right side by E n . Representing the solution 
corresponding to this n-th Landau level by a superscript n, we can then write for the 
positive energy solutions, 

f [ :\y) = ( J % l(e) ) , /<">(*) = ( In ° {0 ) • (22) 

For n = 0, the solution / + does not exist. We will consistently incorporate this fact by 
defining 

/_!(</) = 0, (23) 
in addition to the definition of I n in Eq. (1181) for non-negative integers n. 
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The solutions in Eq. (1221) determine the upper components of the spinors through Eq. 
( TTTj) . The lower components, denoted by x earlier, can be solved using Eq. ([8]), and finally 
the positive energy solutions of the Dirac equation can be written as 



where X M denotes the space-time coordinate. And U s are given by [6j [8] 

/ 4-i (0 \ / o \ 



(24) 



U+{y,n, 







Pz 



E n + m 



( n-l 



(0 



V2neB 



V E n + 



rn 



V2neB 
E n + m 

Pz 



V E n + m 



• (25) 



For the case of positrons which are positively charged negative energy solutions of the 
Dirac equation we have to put Q — — 1 and in this case also we can write their energy as 
in Eq. ( 12~TT) but this time Eq. ( flTl) predicts that the n = solution must only have the 
s = 1 component. Although the dispersion relation of the electrons and positrons become 
different in presence of a magnetic field but they can be written in a unique form as given 
in Eq. fl2T]) . the difference shows up in the spin of the zeroth Landau level state. A similar 
procedure, as used for solving for the positive energy spinors, can be adopted to solve for 
the negative energy spinors and the solutions are: 



(26) 



where 



Pz 



E„ + m 



V 



V2neB 
E n + m 





U0 



( V2ne& ~ \ 

1 n-l{C,) 



I 



E n + m' 

Pz 



E n + m 




V 



(27) 



where £ is obtained from £ by changing the sign of the p^-term. These solutions are 
eigenstates of U x and U z but not of U y . As n^, and Ily do not commute we cannot have 
simultaneous eigenstates of both. 

The solutions of the Dirac equation in presence of a magnetic field are exact solutions 
and not perturbative excitations around the free Dirac equation solutions, which is evident 
from Eq. (1161) . Consequently we cannot put B — > in the final solutions, in Eq. (I25p and 
Eq. f[2Tl) . and expect we will get back the free Dirac solutions. Mathematically in the 
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zero field limit the quantization condition in Eq. (|T6|) fails and in that limit the solutions 
of Eq. (fT5l) becomes indeterminate. Physically we can say that the solutions in Eq. (!25l) 
and Eq. ([271) are specific to a gauge, giving rise to a magnetic field along the z direction, 
and we can at best gauge transform these solutions to obtain equivalent solutions in a 
background magnetic field. The choice of the background gauge does not permit us to 
obtain the free solutions in any limit as the free solutions belong to another gauge orbit, 
namely the pure gauge solutions. 

2.2 The lowest Landau level solutions 

It is previously stated that the n = solution is non-degenerate and in this state we 
have only one solution available for the positive energy and one for the negative energy. 
They are the s — — 1 for the positive energy state and s = 1 for the negative energy 
state, which is evident from Eq. (1251) and Eq. ( [271) . Only in the n = state the wave 
functions are eigenstates of E^, where H z = 27172, and for all other higher Landau states 
the solutions do not have any definite S 2 eigenvalue. In actual calculations when the 
strength of the magnetic field is high we require to work with the n = solutions. We can 
roughly estimate the magnitude of the magnetic field suitable for the n = approximation. 
Suppose we know the typical electron energy in a system is E and the magnitude of the 
magnetic field is B from experimental observations. If it happens that 2neB > E 2 — m 2 
for any positive value of n then from the dispersion relation in Eq. (f2~Tj) we see that p 2 
has to be negative, which is impossible. Consequently when ever 2eB is greater than 
the square of the typical electron energy of the system minus the rest mass square of the 
electron then we have only the n = level contributing to the energy levels and only those 
corresponding wave functions must be used in calculating the other details of the system. 
As an example if the typical electron energy of the system is of the order of lMeV then 
for magnetic field magnitude greater than 10 14 Gauss we must only have the n = level 
contributions in the energy. For lower magnitude of the magnetic field the other Landau 
levels will start to contribute in the electron energy. For a fixed energy of the electron 
and for very low magnetic field magnitude we will have many possible Landau levels. 

2.3 Ortho-normality of the spinors and their completeness 

Using the relation 



where 5 n ^ m = 1 when n = m and zero otherwise and a is dimensionless we can calculate 
the ortho-normality of the spinors. The ortho-normality of the spinors in the present case 
has to be modified as the spinors have explicit co-ordinate dependencies. Using Eq. ( 1281) 
it can be shown in a straight forward fashion that, 



dyU}(y,n,i\ t )U s/ (y,m,i^) = / dyV*(y,n,p i( )V s >(y,m,p k ) = 5 nim 5 8i8 > " ,(29) 




(28) 



2E, 
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and 

/oo roc 
dyUj(y,n,^ t )V sl (y,m,-p k ) = dyV](y,n,-p k )U s ,(y,m,p y ) = 0. (30) 
-oo J — oo 

Except the integration over y and the appearance of the Landau levels the above relations 
closely resemble the corresponding relations in free-space. The above relations fix the 
normalization of the spinors. We will rederive the normalization constants of the spinors 
when we quantize the theory in section HI 

Using now the solutions for the U and the V spinors from Eqs. (12"oT) and (12"T|) . it is 
straight forward to verify that, 

/OO j . 

dy ]T (U s (y, n, n) N U^(y, ™, ft ) + V s N (y, n, -p^)V s N \y, m, = 6 n , m l (31) 

Here (y,n,j^) and V s N (y,n, — p^) are the normalized spinors and 1 is the unit 4x4 
matrix. If the Landau levels of the two spinors are the same then we have, 

]T (U s (y, n, p^UKy*, n, ft,) + V s (y 1 n 1 -p k )V}(y^,n,-p i ) 



( v 2 + 2neB \ r 

= 1 + % ^ xdia § Vi(e)Vife), / n (£)/n(£*),4-i(O4-i(£0, 4(04(6 

(32) 

where 'diag' indicates a diagonal matrix with the specified entries, and £ and involve 
the same value of p x . A sum over the Landau levels for spinors situated at different y 
co-ordinates gives, 

OO 

Y / T,{ U ^y^^n)U^(y,,n,p y ) + V s N (y,n,-p il )V s N \y^n,~p y )) = 5{y - Vi ) 1, (33) 

n=0 B 

where we have used the result of Eq. (1201) . The two equations in Eq. f[3"Tj) and Eq. fl3"31 
stands for the completeness relations for the spinors in the present case. 

3 Spin-sum of the Dirac solutions in an uniform back- 
ground magnetic field 

In this section we derive the spin-sum J2 S U s (y, n, j\,)U 8 (y*, of the solutions of the 
Dirac equation in presence of a magnetic field. The two spinors in the above sum can have 
two different position coordinates in general and so their spatial dependence is explicitly 
shown to be different. From the nature of the solutions as given in Eq. ( 125]) we see that 
T, s U s (y,n,j\)U s (y*, n,j\) can be written as: 

_ 1 n 

Pu(y,y*,n,j\) = 2_,U s {y,n,p k )U s {y itl n,p u ) = — — ■ £(0-0(60 T i,j ( 34 ) 
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The spin-sum of the product of the spinors, J2 S U s (y,n,p jt )U s (y i ,,n,p y ) will give rise to a 
4x4 matrix whose elements will be contain where i, j runs from n — 1, n. 

If these terms as }(£*) are taken as common factors then the whole 4x4 spin-sum 
matrix can be represented as a sum of terms containing the products of times 
the corresponding 4x4 matrices called T^j. 

Using the dispersion relation E 2 = pi + m 2 + 2neB, T n , n can be written as [7] , 



/ 

(E n + m) 



V o - Pz o 







Pz 



<E n -m) j 



In the 2x2 notation the above matrix can be written as, 



T 



E u 



+ m 



1(1-^ 











5(1-0-3) \ 



+ Pz 



■|(l-0 3 ) 



fU - 03J 








1(1 -03) 



(35) 



(36) 



where 03 is the third Pauli matrix. In the 4x4 notation Eq. (|36l) can be written as, 



T 



-[m(l - E.) + £ n ( 7 u + 7 V) - p 2 ( 7 V + 7 d )] , 



;H1-S 2 )+^|| +A75], 



(37) 



where <r 2 = ij lr y 2 . In the last equation f\\ = p°j + p 3 7 3 and ftu 
75 = «7°7 1 7 2 7 3 . In our case po = E n . 

In a similar way T n _i n _i can be written as: 



T 

- 1 n— l,n- 



In the 2x2 notation the above equation looks like 

'5(1 + 03) N 



5(1 + 03) 



P°73 + P 3 7o and 



/ (£„ + m) 


-p z 


M 













-{E n -m) 





V 





/ 



•Ml- 1,71- 



+ m 



-5(1 + 03) ; 
^ 

i(i+ as )^- 

In the 4x4 notation the above equation becomes 

1 



+ Pz \ 1 



Kl + 03) 



5 1 + 03 



T 

- t n— 1,71- 



m(l + £ 2 ) + £ n ( 7 ° - 7 5 7 3 ) + p 2 ( 7 5 7° - 7 3 )] 



m 1 



l + £ 2 )+^l -^75]. 
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(38) 



(39) 



(40) 



From the matrix multiplication in the left hand side of Eq. ([Mil it can be seen that 
T n _i n is given as, 



f 








1 \ 

















-1 








V o 








o J 



T n ^ 1>n = V2neB 

In the 2x2 notation the above equation looks like 





(41) 



T n _!, n = V2neB 



(cji + ia 2 ) 



2(0-1 + 101) 




(42) 



Here o\ and o-i are the first two Pauli matrices. When converted back to the 4x4 notation 
the above equation becomes, 



T n -i, n = -^V2nei3(7i + ij 2 ) 



(43) 



Similarly T„ in _i is given by, 



T n ^i = V2neB 
In the 2x2 notation the above equation looks like, 



( ^ 

10 



V -1 J 



(44) 



V -3(^1 - ia v 





which when converted back to the 4x4 notation becomes, 

172) • 



T n -i, n = -^V2nei3(7i - 



(45) 



(46) 



Supplying the values of T itj s from Eq. (j37jh Eq. (j40T), Eq. (J43]) and Eq. (J46]) to Eq. ([34 
we get the result: 



Pu(y,y*,n,; 

where, 

Su(y,y*,n,p !t ) -■ 



= ^U s {y,n,p H )U s {y^n,p H ) 



[E n + m) 



Su(y,y*,n,j\) , (47) 



[mil + £,) -^ |75 }/ n _i(e)/„-i(e*) 
{m(l-S z )+^+^ ||75 }j n (e)/n(^) 

\/2^B( 7l - i l2 )I n {£)I n _ x {Q - v / 2^B( 7l + i7a)/„_i(0/n(C* 



(45 
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Similarly, the spin sum for the V^-spinors can also be calculated, and we obtain: 

Pv(y,V*,n,H) = ^V a (y,n,Q t )V a (y,n,R,) = rS v {y,y in n,j\) , (49) 

where, 

1 



Sv(y,V*,n,i\) 



{-m(l + S 2 ) + ^ - ^ ||75 } 



+ V2neS( 7l - i 72 )7 n (0/„-i(6) + V2ne/3( 7l + 172) I n _i (£)/„(& 



(50) 



One important property of the above spin-sums is that, 

Puiv.y^n,^) = -P v {y,y^n, , (51) 
which is similar to the result in vacuum. 

4 Quantization of the electron field in presence of an 
uniform background magnetic field and the elec- 
tron propagator 

In this section we will use the spin-sum results in writing the electron propagator in 
presence of an external uniform magnetic field. But before doing so we will first write 
down the QED Lagrangian for the electron in presence of a background magnetic field. 

In presence of a background magnetic field we can decompose the photon field as 
follows: 

A»(x) = A£(x) + A£(x) , (52) 

where A^x) is the dynamical photon field which will be quantized and A B (x) is the 
classical background field which gives rise to the magnetic field. If the uniform background 
classical magnetic field is called S then we must have: 

B = Vx A B (x) , (53) 

where A B (x) = (0, A B (x)). In presence of the background magnetic field we can also 
write the field strength tensor as: 

F^(x) = F^(x) + F^, (54) 

where F£ v {x) = d^A^(x) — d v A^(x) and F$ = d l A B (x) — d j A l B (x) is a constant as given 
in Eq. §3J. 
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The QED Lagrangian can be written as: 

£ = $(z 7/1 £>" - m)V> - \f^F^ , (55) 

where = — ieA^ is the covariant derivative of the fermion fields. The QED La- 
grangian can also be written as: 

C = i) [ 7/ ,rP -m]ip + e^MS - \f^F^ , (56) 

where IP = id^ + eA^ is the kinetic momentum of the fermions in presence of the back- 
ground field. The first term of the Lagrangian contains no dynamical photon dependence 
but it depends upon the background magnetic field through II and this part of the La- 
grangian gives rise to the Hamiltonian of the electron in presence of the magnetic field 
used in Eq. (jlj). The equation of motion which we obtain from the first term of the above 
Lagrangian is in fact the Dirac equation in presence of a magnetic field which we solved in 
section [21 Consequently the most important effect of the background magnetic field is to 
modify the solutions of the Dirac equation. The interaction term of electrons and photons 
remains the same as in normal QED. The free fermionic part of the Lagrangian in Eq. (1561) 
is also important for the definition of the propagator of the electron and in the next part 
of this section we will find out the expression of the electron propagator in presence of 
the background magnetic field. Before we calculate the electron propagator we quantize 
the theory. The photons do not interact with the magnetic field and consequently their 
quantization procedure is the same as in normal QED. 



4.1 Quantization of the electron field 



Since we have found the solutions to the Dirac equation, we can now use them to construct 
the fermion field operator in the second quantized version. For this, we write 



*{X) 



EE 

s=± n=0 



dPxdPz 
2ttD 



/ s (n,^)e-^C/ s ( Z /,ri, H ) + ^(n,^)e i ^y s (y,n, H )l . (57) 



Here, / s (n,]^) is the annihilation operator for the fermion, and /j(7i,J\) is the creation 
operator for the antifermion in the n-th Landau level with given values of p x and p z . It is 
to be noted that the wave-functions of the electron used in Eq. (157)) are not free-particle 
solutions and they never tend to the free-particle solutions in any limit. As the the 
positive and negative frequency parts of the solutions are as free-particles so the notion 
of a particle and anti-particle is unambiguous in the present circumstance. The creation 
and annihilation operators satisfy the anti-commutation relations 



$ss'$nn'd(Px -p' x )5(j) z - p' z 



(5f 



and a similar one with the operators / and f>, all other anti-commutators being zero. The 
quantity D appearing in Eq. (157)) depends on the normalization of the spinor solutions, and 
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in this section we will rederive the normalization of the spinors calculated in subsection 
12.31 The factor of 2n multiplying D helps to associate D with the normalization constant 
found in subsection 12.31 Once we have chosen the spinor normalization, the factor D 
appearing in Eq. (157)) is however fixed, and it can be determined from the equal time 
anti-commutation relation 



5\x - x,) . 



(59) 



Plugging in the expression given in Eq. (1571) to the left side of this equation and using the 
anti-commutation relations of Eq. (15H1) . we obtain 

dp x dp z 



( e -iM*-s* ) e - ip ^-^U s (y,n, n )Ut(y„n,p k ) 
+e ^-^e i ^-^V s (y,n, R )V](y^n, n ) 



(60) 



Changing the signs of the dummy integration variables p x and p z in the second term, we 
can rewrite it as 

dp x dp z 



^(X),^pQ 



EE 



^^-^e-^-^^U^y^^^Ky^n,^) 



+V s (y, n, -fl,)V?(y*, n, . 



(61) 



At this stage, we can perform the sum over n in Eq. (161)) using Eq. (132)) and Eq. ( 1201) 
which gives the (^-function of the ^-coordinate and perform the integrations over p x and 
p z to recover the ^-functions for the other two coordinates as well, provided 

E n + m (2vrD) 2 (2ir) 2 ' 1 ' 

In this way we get back the same value of the normalization of the spinors which we 
obtained in subsection 12.31 Putting the solution for D, we can rewrite Eq. ( 1571) as 



EE 

s=± n=0 



dp x dp z E n + m 



2tt 



2E n 



(63) 



x [/ s (n,^)e-^C/ s (y,r i , H ) + ^(n, H )e^K(y,r i ,^y . 

The one-fermion states are defined as 

\n,p,,s) = 0/1(71,^)10) . (64) 

The normalization constant C is determined by the condition that the one-particle states 
should be orthonormal. For this, we need to define the theory in a finite but large region 
whose dimensions are L x , L y and L z along the three spatial axes. This gives 

V L X L Z 

Next we calculate the electron propagator in presence of an uniform background magnetic 
field. 
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4.2 The electron propagator 

The electron propagator is given by, 

iSfiX-XJap = (0\T{ij;(X) a lp(XM\0) 

= 6(t - Q(Q\^(X) a ^(X^\Q) - 9{U - t)(0|^)^(X) Q |0) , 

(66) 

where T is the time-ordered product and 0(A) is the step-function which is unity when 
A > and zero other wise. The step function can be represented in an integral form as: 

, (67) 

-co 2tt uj — le 

where e is an infinitesimal parameter. Using Eq. ( |63i) we can write, 

(OmX)^(XM\0) = ffj^(^^)e-^- x ^ 

x U 3ta (y, n, p^U s ,p(y*i n > Py) > 

(68) 

and using Eq. (HTj) and suppressing the spinor indices the above equation can also be 
written as, 

W(XMXOIO) = tj d ~j0 (^f) e-^-^Puiy^n^) . (69) 
In a similar way it can be shown that, 

mX^(X)\0) = £ / d j0- e^- x ^P v (y,y^) , (70) 



where Py is given in Eq. (149]) . Using the above results in Eq. (I66I) and utilizing the form 
of the ^-function in Eq. (|67|) we can write, 



CO „ 

iS*(X - XJ = i^, J 



dpx dp z duj ( E n + m 



n=0 

" g— icj(t— U) —ip- (X^-X^ ) iu>(t—U)+ip- (X^ ~X^ ) 

x : P u (y,y in n,p H ) : P v (y,y^n,i\) 



(2tt) 3 V 2E,, 



u — le u — le 



.(71) 



Changing the signs of p x and p z in the second term of the integrand and using Eq. (joT 

we get, 



(72) 



n=0 

X 



i(w+B„)(t-t*) e i{u+E n )(t-ti,) 

: + : 

uj — le uj — te 
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Now appropriately doing the integration over to we get, 



oc 



is§(x-x+) = iJ2 I dPxdPzdp ° c - ip - (x ^ XA) s u(y,y*, n ,i\) 



n=0 
oo 



(27r) 3 Pq — p 2 — m 2 — 2neB — ie ' 



\- / dp x dp z dp . (x ^_ Xa) S u {y,y^n,p H ) 
^ I (2tt) 3 p\ - m 2 - 2neB 



i 

n=0 
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where Sjj(y, y±, ^,f\,) is given by Eq. ( 1481) and p?< = pi — pi- It is to be noted that the pole 
of the propagator is now dependent on the Landau levels as it should be in an uniform 
background magnetic field. The form of the propagator suggests that it is not translation 
invariant and so it cannot be written down completely in Fourier space. 



5 A discussion on gauge dependence 

As we have solved the Dirac equation in presence of a uniform background magnetic field 
using a particular gauge, as given in Eq. (TjQ) with b = 0, the solutions are dependent on the 
gauge choice. The spinor solutions are themselves not physical observables and so they can 
be gauge dependent. But not all the results discussed in this article are gauge dependent. 
The energy of the electron as given in Eq. (fT7|) is not a gauge dependent quantity, any 
gauge we choose we will get the same dispersion relation of the electrons. The special forms 
of the ortho-normality relations as given in section 12.31 are gauge dependent as the results 
contain the functions which has p x which is not a gauge invariant quantity. The spin-sum 
also depends on the particular gauge we work in and the above results will be different if we 
had chosen another gauge to represent the magnetic field. But actual calculations yielding 
physical quantities like scattering cross-section or decay rates must be independent of the 
choice of the background gauge field. We can see the gauge invariance of the physical 
quantities in a heuristic way. If we had chosen the gauge specified in Eq. §2§ with c = 
instead of the gauge in Eq. (TjQ) with b = then the solutions of the Dirac equations as 
specified in Eq. (125]) and Eq. (1271) should have been the same except all the y should be 
replaced by x and p x should be replaced by p y inside the spinors and the free-particle 
part should contain e tp ' x * instead of e tp ' X)i . A similar replacement should yield the new 
spin-sums and the propagator. Consequently the quantities calculated in these two gauges 
differ by the way we name the x and y coordinate axes. But in calculations of scattering 
cross-sections and decay-rates we always have integrations over x,y, z coordinates at each 
vertex and consequently the end results will not depend upon which gauge we started 
with. 

The above discussion highlights the fact that most of the quantities calculated in this 
article using the exact solutions in presence of the magnetic field rely heavily on our choice 
of the vector potential. All the solutions of the Dirac equation in presence of a uniform 
magnetic field along the z direction obtained by using various vector potentials will be 
different but are related by gauge transformations. It is to be noted that the free Dirac 
solutions can also be gauge rotated where the gauge fields are pure gauge configurations. 
As there is no connection between the gauge configurations giving rise to a magnetic field 
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along the z direction and pure gauge fields so we do not get back the free Dirac solutions 
limit of the exact solutions in a magnetic field. 



6 Conclusion 

In this article we solved the Dirac equation in presence of a background uniform magnetic 
field specified by a particular vector potential. The dispersion relation of the electron is 
seen to change from its form in the vacuum and we see the emergence of Landau levels 
designating the quantized nature of the transverse motion of the electrons. The solutions 
of the Dirac equation are dependent on the Landau levels, the energy of the electron is 
seen to be degenerate except the lowest Landau level energy. It is seen that there is no way 
to get back the free Dirac solutions from the exact solutions in presence of the magnetic 
field by letting the field strength to go to zero in the solutions, a fact which is related to 
the gauge invariance of the system. Using the appropriate spinors in a magnetic field the 
ortho-normality and completeness of the spinors were worked out in section 12.31 which 
closely resembles the corresponding results in vacuum. The spin-sum of the solutions are 
derived explicitly using the exact solutions of the Dirac equation in a magnetic field. The 
theory is quantized and with the quantum field operators the propagator of the electron 
in presence of an uniform background magnetic field is calculated. Some thing similar to 
our derivation of the electron propagator was presented in [9] where the authors worked 
in the chiral representation of the Dirac gamma matrices. But the presentation of the 
expression of the propagator was not compact and nor the authors in [9] calculate the 
spin-sum explicitly. As most of the quantities calculated in this article depend on the 
choice of the vector potential giving rise to the magnetic field so the gauge invariance of 
the calculations become less transparent. In the penultimate section we discuss about the 
gauge invariance of the calculations in presence of a magnetic field and show that although 
the spin-sums and propagators may not be gauge invariant but physical quantities like 
scattering-cross sections and decay rates of elementary particles in presence of a magnetic 
field can be gauge invariant. 
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